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Abstract 



The temperature dependence of conductance resonances is used to mea- 
sure the evolution with the magnetic field of the average level spacing Ae of 
a droplet containing ~ 30 electrons created by lateral confinement of a two- 
dimensional electron gas in GaAs. Ae becomes very small (< 30/ueV) near 
two critical magnetic fields at which the symmetry of the droplet changes 
and these decreases of Ae are predicted by Hartree-Fock (HF) for charge 
excitations. Between the two critical fields, however, the largest measured 
Ae = 100/ieV is an order of magnitude smaller than predicted by HF but 
comparable to the Zeeman splitting at this field, which suggests that the spin 

degrees of freedom are important. 
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In a recent Letter, Klein et al. (l[ showed that a droplet of electrons created by lateral 
confinement of a two-dimensional electron gas undergoes changes in its symmetry at two 
critical magnetic fields. On one hand, the exchange and the confinement energies favor a 
compact electron distribution, on the other hand, the Coulomb repulsion favors a diffuse 
occupation. The magnetic field B alters the balance between these two effects, causing 
the ground state (GS) to change ||. Correlations are also expected to play an important 
role ||. Such symmetry breaking is expected to be accompanied by low energy excitations. 
In particular, Hartree-Fock (HF) [fjj calculations, which predict these changes in symmetry, 
also predict that these are accompanied by a decrease of the average excitation energy of the 
droplet. We report, here, measurements of the average level spacing Ae(B) which confirm 
these qualitative predictions. Between the two critical magnetic fields, however, the largest 
measured Ae is much smaller than predicted by current theories for charge excitations. The 
maximum Ae(B) is close to the Zeeman energy, which suggests that spin excitations are 
important. 

The device we use consists of a heavily n-doped GaAs substrate covered by a first layer 
of AlGaAs and a top layer of GaAs ||, both layers grown by molecular-beam epitaxy. The 
strong electric field created by the the band offset between AlGaAs and GaAs, as well as 
the positive voltage applied to the substrate on the bottom, creates a two-dimensional 
electron gas (2DEG) at the AlGaAs/GaAs interface. Ti-Au electrodes are deposited on 
the top surface lOOnm above the 2DEG and fashioned by electron beam lithography. A 
negative voltage maintained constant during the experiment, is applied to the top electrodes. 
This confines the electrons laterally in a potential that is approximately a two-dimensional 
harmonic oscillator characterised by the energy huo [§J. In addition, the electrodes define 
two potential barriers, through which electrons must tunnel to enter and exit the droplet. 
The 2DEG regions outside the left and right barriers form the leads. The current through 
the droplet is measured as a function of the voltage V&, as well as that between the left and 
right leads V\ r . 

Figure 1(a) illustrates the conductance G through the droplet for VJ r < ksT/e, where G 
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is independent of V\ r . Sharp peaks arise from resonant tunneling: the iVth peak occurs at 
a voltage such that eVj^ is proportional to (E N — E N _i), the energy difference between 
a (iV)-particle and (N — l)-particle GS. At high temperatures, this difference is dominated 
by the Coulomb charging energy of the droplet. At very low temperatures, however, there 
is a quantum mechanical contribution to (E^ — -Ejv-i) arising from the confinement of the 
droplet to a small region of space. When the GS of the droplet changes, (En — -Ejv-i) 
changes as well and thus Vjv provides a spectroscopic probe of the GS energy of the droplet. 

Figure 1(b) illustrates how this addition spectroscopy is done by following the 

-B-dependence of the peak position at the base temperature of our dilution refrigerator 
corresponding to an electron temperature of 50mK (see caption of Fig. 2). At high B, the 
kinetic energy of the electrons is quantized in Landau levels (LLs) of index n = 0,1, .... The 
large structure in V^(B) below 1.5T corresponds to the transfer of electrons from higher to 
lower n LLs. The change in behavior at 1.5T indicates a different phenomenon. Above this 
field all the electrons are in the lowest LL (i.e. all electrons have the same kinetic energy) 
and the jumps in Vjv are caused by the flips of electron spins. The spins flip, not because 
of their magnetic moment, which is small in GaAs (g = —0.4), but because of the Coulomb 
interaction 0,|TJ. The latter scales as e 2 /e£e where £b — \Jfrc/eB is the magnetic length 
and e is the dielectric constant. Increasing B decreases is, which increases the Coulomb 
repulsion, causing the electrons to spread out |j,fE|. 

Klein et al. |IJ have shown that just below B c (indicated in Fig. 1(b)) the GS of the droplet 
is a singlet. Their experiment shows that for a droplet with a small number of electrons 
iV ~ 30 and a large hu> ~ 2meV, the GS has a compact charge distribution at low fields. 
The compact occupation of the lowest LL is achieved when the charge density corresponds 
to filling fraction (the ratio of electron to flux-quantum density) equal to 2 throughout the 
droplet, i.e. when both spin states are equally occupied. Above B c , the growth in the 
Coulomb repulsion causes the charge to spread out, causing the electrons to flip their spins 
in order to minimize the loss of exchange energy. Thus the transition at B c results from a 
change in symmetry from zero to finite total spin. 

3 



A different kind of symmetry breaking occurs at high fields. Above 1.9B C all spins are 
polarized and the droplet must find a new way to spread out as B increases. In HF, the GS 
between 1.9B C and B r = 2.5B C (Fig. 1(c)) is the so-called maximum density droplet (MDD) 
|7]j2|], the compact charge distribution for a spin polarized droplet with filling fraction 1 
throughout the droplet. In HF, the symmetry change at B r is analogous to a liquid- gas phase 
transition and the order parameter is the charge density. Klein et al. observe experimentally 
that B r is lower than predicted by HF, as seen by comparing Fig. 1(b) and (c), suggesting 
that correlations play an important role in this higher field symmetry breaking H,[|. 

HF predicts a dramatic decrease in the average excitation energy Ae of the droplet near 
B c and B r . One way we determine Ae is from the cross-over from single to multiple level 
transport as T is increased 0. We carefully select conductance peaks that have an expo- 
nential tail at base temperature (Fig. 1(a)). This indicates that the peak shape is dominated 
by the thermal broadening of the energy distribution of the electrons in the reservoir || and 
that quantum fluctuation effects |K| are comparatively small. In contrast, a Lorentzian tail 
is the characteristic signature of the regime where the latter effects are important. Also, we 
measure the temperature dependence at B fields well separated from the steps in Fig. 1(b), 
for at the cusps associated with these we expect degenerate GS's. 

When ksT < Ae, the current is limited by a single quantum level and the conductance 
peak profile is given by the derivative of the Fermi-Dirac distribution function. The data in 
Fig. 1(a) are well fit by the formula [11J 

Tat is the tunneling matrix element, and the factor a converts a change in V& to a shift 
in the electrostatic potential of the droplet RI2| . Thus, the amplitude of the peak G max 



decreases as 1/T with increasing T. However, when /c#T becomes larger than Ae, excited 
states as well as the GS participate in the conductance. The number of levels participating 
then grows as T/Ae, but each channel still contributes a weight that varies as 1/T, so the 
total conductance becomes temperature independent. The cross-over of G max from 1/T 



to constant provides a measure of Ae. Beenaker [|IT] has calculated the line shape in the 
multi-level regime assuming a constant Ae and a constant for all levels in the excitation 
spectrum of the (iV)-electron droplet. He finds that the following is a good approximation: 

The cross-over temperature of G max from 1/T to constant is at Ae/2. 

Concomitant with this cross-over of the peak amplitude, the peak profile also changes 
subtly, and this shows up in measurements of the full width at half maximum (FWHM) 
||. The FWHM is proportional to T for both single- and multiple-level transport, but the 
FWHM is larger in the multiple-level regime. The increase of the FWHM also provides an 
indication of when fcgT is comparable to Ae. 

Figure 2(a) shows the T-dependence of the conductance peaks at B/B c = 1.6. The 
inverse of G max as a function T is plotted in Fig.2(b). The cross-over from 1/T to constant 
is quite clear. The FWHM as a function of T for one of the peaks in Fig. 2 (a) is plotted in 
Fig. 2(c). At T ~ 0.5K where the amplitude cross-over takes place, the FWHM also deviates 
from its low T-linear behavior. 

Small differences from the model of Eqs.(Q) and (Q) are observed in Fig.2(b). We observe 
that 1/G m ax decreases slightly with increasing T above 600mK. This can be explained by 
taking into account the variation of T and Ae with excitation energy for fixed N. One 
expects that levels of increasing excitation energy have smaller Ae and larger T because of 
deviations from parabolicity of the confining potential and because of narrowing of the two 
tunnel barriers at higher energy. 

Figures 2(d-f) show analogous measurements at B/B c = 2.1 close to B r . Clearly Ae is 
smaller near the field at which the symmetry breaking occurs. The T-dependence of the 
FWHM is also consistent with a small Ae. The FWHM is linear in T down to the lowest 
temperature measured, with a slope larger than the one measured in Fig.2(c), suggesting 
that the droplet is already in the multiple-level regime at the lowest temperatures. 

Figure 3(a) shows the 5-field dependence of Ae extracted from data like those in Fig. 2 



for devices with 30 and 50 electrons in the droplet. Although the values of B c differ by 0.2T 
between the two cases, the experimental results for the two droplets are consistent when 
plotted as a function of B/B c . The level spacing nearly vanishes at B/ B c = 1 and 2.2 where 
the two symmetry changes occur. Between these two critical fields, Ae reaches a maximum 
value of 100/xeV. 

The inset of Fig. 3(b) illustrates the HF prediction for the level spacing. We have calcu- 
lated the gap between the GS and the lowest excited state as a function of B. The latter 
energy vanishes at level crossings, which occur near steps in Fig. 1(c), and rises to a maxi- 
mum between steps. We have plotted, in the inset of Fig.3(b), the lines that join all these 
local maxima. In the experiment, we have also measured Ae between steps. The overall 
shape of Fig. 3 (a) is reproduced by HF: Ae vanishes near B c , rises to a maximum, and then 
falls again precipitously near B r . However, the largest values of Ae predicted by HF are ten 
times larger than observed experimentally. 

Another way to measure Ae is to use tunneling excitation spectroscopy (TES) p!3|JT^] . In 



this scheme, the differential conductance is measured as a function of V\ r and the quantum 
levels are seen as peaks. Two spectra are shown in Fig.4. At B / B c = 1.6, where Ae is largest, 
the level spacing measured this way agrees well with that measured from the T-dependence. 
However at B/B c = 2.1, close to B r , TES gives a value somewhat larger than does the 
temperature dependence, although the value is much smaller than at B/B c = 1.6. TES 



measures the single-particle excitations [|15| while the T-dependence also includes many- 
body excitations. This may be the reason for this discrepancy. 

The observation of very small Ae near B c and B r confirms the interpretation of Klein et 
al. [[H that symmetry breaking occurs at these critical fields. However the observation that 
the largest Ae is much smaller than the value predicted by HF calls into question the nature 
of the low-lying excitations of the MDD. 

The existing models that describe the low-lying excitations near the MDD can be sep- 
arated in two categories depending on whether they consider charge or spin excitations. 
Consider, first, the former, which assumes that both GS and excited states are spin polar- 



ized. This class includes the HF calculation of Chamon and Wen J2|, which uses as a basis 
the single-particle states of the symmetric gauge, without level mixing. The single-particle 
states are labeled by an angular momentum index m > 0, and they represent circular 



orbits of radius £sy2(m + 1). In the MDD GS, all the innermost orbitals are occupied 
(m = 0, 1, ..,N — 1) and the charge density has the spatial distribution corresponding to 
filling fraction 1 uniformly over the droplet. This compact charge distribution is incompress- 
ible. As B increases, however, the HF charge distribution approaches the classical dome-like 
shape corresponding to a state which is compressible throughout the droplet. The first step 
in this transformation is the edge reconstruction at B r where holes are first introduced in 
the interior of the droplet. This transition, from an incompressible state to a compressible 
one, gives rise to the abrupt decrease of the level spacing at B r . The single-particle level 
spacing in the compressible state can be evaluated from the expression 2/i 2 /m*r 2 g,, ~ 40/ieV 
where m* is the effective mass, g s = 1 is the spin-degeneracy and r = 200nm is the radius 



of the MDD. A better approximation than HF is the model of Oaknin et al. (16| in which 
charge magnetoexcitons are the low-lying excitations. These are single electron-hole pair 
excitations of the MDD which correspond to moving an electron from the interior to the 
exterior of the droplet. While such excitons give rise to excited states in the MDD, they 
become stable in the GS above B r . Above B r , as in HF the electron occupancy is reduced 
within a few £g inside the edge of the MDD. Other models that incorporate correlations 
include the work of Kamilla and Jain fl7|,|3|l that study excitations of non-interacting com- 



posite fermions. In all the models we have discussed, the level spacing in the MDD state 
is of order e 2 /e£s, which is approximately 7meV near 3T, the 5-field at which the MDD is 
formed for the data in Fig.3(a), and is therefore much larger than the Ae we observe. 

In contrast, the Zeeman energy is 75/ieV at 3T (using g = —0.4), a value closer to our 
measured Ae. Several calculations have appeared recently on spin-wave like excitations [18| 



in a droplet ||19|| . By canting the spins of the electrons gradually over the MDD, the droplet 
reduces the cost in exchange energy of having two neighbouring electrons with opposite spins. 
Such states may also be a better description of the GS near B r . Even if the MDD is a good 



description of the GS, excited states that involve spin excitations may have lower energy 
than those involving charge excitations, which cost both confinement and exchange energies. 
In particular, a uniform rotation of all spins costs no exchange or confinement energy at all. 
The total spin quantum numbers of the spin-polarized MDD are S 2 = {N/2 + l)N/2 and 
S z = N/2. The first excited state has the same S 2 as the MDD but S z = N/2 — 1 and the 
energy gap is then given by the Zeeman energy. This excited state can also be obtained 
by including vertex corrections in the Coulomb interactions [20]. While, this gives Ae of 
the right order of magnitude, the actual excited states may be more complex, involving an 
admixture of both spin and charge excitations depending on N and the Zeeman energy. 

We have plotted in Fig.3(b) the lowest Ae(£>) predicted for either spin excitations or 
charge excitations. The upper branch of the trapezoidal shape in Fig. 3(b) is the Zeeman 
energy. The two abrupt decreases at B c and B r are the low-lying charge excitations predicted 
by HF near the fields at which symmetry changes. 

Our measurements confirm the predictions of HF that Ae vanishes near two critical fields 
at which changes in the symmetry of the droplet occur. Between these two critical fields, 
however, the spacing is of the order of the Zeeman energy suggesting that the low-lying 
excitations are those of the total spin of the droplet rather than charge excitations included 
in HF and some other models. 
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FIGURES 

FIG. 1. (a) Linear conductance through the droplet as a function of the voltage on the bottom 
gate Vb at B = 2.7T. (b) Variation with B of the position of the peak near 0.160V. (The -B-sweep 
rate is 0.025 G/sec.) We count 13 steps between B c and 2B C which means that N = 26 or 27. (c) 
HF calculation of the resonant energy as a function of B for a droplet containing N = 27 electrons 
and Twoq = 2.1meV. The model has the same number of steps and the same B c as the experiment. 
Note that the height of the steps is larger in HF (c) than in the experiment (b) suggesting that 
the excitation gap is also larger in HF. 

FIG. 2. Conductance vs. Vj, for T ranging from lOOmK to 600mK in increments of lOOmK, 
(a) measured at B/B c = 1.6 and (d) at B/B c = 2.1. (b) and (e) Inverse of the conductance peak 
vs. T for the peaks in (a) and (d) respectively. The cross-over from 1/T to constant determines 
Ae/2. (c) and (f) Full width at half maximum vs. T for one peak in (a) and (d) respectively. The 
error bars are determined by comparing the behavior of other peaks. In (c), the deviation from 
the low T straight line gives a similar value of Ae/2 as in (b). The width and amplitude at base 
temperature indicate that the 2DEG is at 50mK, whereas the base temperature of the dilution 
refrigerator is 25mK, measured by nuclear orientation thermometry. 

FIG. 3. (a) Ae obtained using the analysis described in the text as a function of B/B c . The 
open circles are for a droplet containing N ~ 30 electrons and the solid circles are for one containing 
N ~ 50 electrons, (b)inset: HF calculation of the S-dependence of the gap between the GS and 
the first excited state. The gap vanishes near each of the steps in Fig. 1(c) and rises to a maximum 
between. The solid line joins all the local maxima between steps. The dashed line is the Zeeman 
energy. Main: independence of the smallest of these two curves shown in the inset. The result 
illustrates Ae(-B) for both charge or spin excitations. 

FIG. 4. Tunneling excitation spectra, dI/dVi r as a function of V/ r at (a) B/B c = 1.6 and (b) 
B/B c =2.\. U ~ 0.5meV is the Coulomb gap and Ae is the level spacing. A factor (5 = 0.7 converts 
the x-axis scale to meV units [12]. 
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